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Abstract 

The isoscalar octupole response of a heavy spherical nucleus is analyzed 
in a semiclassical model based on the linearized Vlasov equation. The oc- 
tupole strength function is evaluated with different degrees of approximation. 
The zero-order fixed-surface response displays a remarkable concentration of 
strength in the Ifiu and 'ifiuj regions, in excellent agreement with the quan- 
tum single-particle response. The collective fixed-surface response reproduces 
both the high- and low-energy octupole rsonances, but not the low-lying 3~ 
collective states, while the moving-surface response function gives a good 
qualitative description of all the main features of the octupole response in 
heavy nuclei. The role of triangular nucleon orbits, that have been related 
to a possible instability of the spherical shape with respect to octupole-type 
deformations, is discussed within this model. It is found that, rather than 
creating instability, the triangular trajectories are the only classical orbits 
contributing to the damping of low-energy octupole excitations. 
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Keywords: Vlasov equation, isoscalar octupole resonances, octupole surface modes 

I. INTRODUCTION 



It is well known that there is an intimate connection between the shell structure in 
quantum systems like nuclei and metallic clusters and the properties of classical trajectories 
within these systems (see e.g. [1], p. 579). In particular, for nuclei it has been argued that 
closed orbits of triangular shape might lead to an instability of the spherical shape against 
octupole-type deformations in the region beyond ^osp]^ j^j^j^ p. 560. Here we would like to 
investigate in detail this possibility by using a semiclassical theory of nuclear response based 
on the linearized Vlasov equation [2,3]. Instabilities are expected to show up as some kind 
of pathological behaviour in response functions (vanishing eigenfrequencies and diverging 
response) and the semiclassical theory of [2] and [3], that has already been shown to give 
good qualitative results for lower multipolarities [4-6], is an ideal tool to study the role 
played by classical trajectories in determining the response of large quantum systems. Of 
course quantum corrections are expected to modify the results of this theory [7,8], especially 
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at very low energy, however having a clear picture of what should be expected already at 
the classical level might help in making progress. We study the isoscalar octupole response 
of a sample "nucleus" of A = 208 nucleons contained in a squarc-wcU potential of radius 
R — ro^a, with tq = 1.2 fm. This is done with different degrees of approximation: in a first 
approximation, discussed in Sect. II, the interaction between nucleons is neglected, while 
in Sect. Ill this interaction is taken into account in a separable approximation. Finally in 
Sect. IV the effect of surface vibrations on the octupole response function is included. 

Because of an interesting scaling property of the zero-order scmiclassical response, the 
^-dependence can be factorized and becomes trivial. Thus, although we perform explicit 
calculations only for A — 208, our results on stability apply also to the region beyond ^°^Pb. 



II. SINGLE-PARTICLE RESPONSE (FIXED SURFACE) 

Following Ref. [1], the connection between shell structure and classical trajectories is 
most easily illustrated by considering the quantum energy levels of a particle in a spherical 
potential (no spin-orbit interaction). In this case the energy levels depend on the radial and 
angular momentum quantum numbers rir and I only. For sufficiently large values of these 
quantum numbers, the energy difference between two neighbouring levels can be approxi- 
mated as 

6(n;, - e{nr, I) ^ « -nr)^ + (/' - 0§ ■ (2-1) 

Now the important point is that the derivatives de/drir and de/dl are essentially the 
derivatives of the classical Hamiltonian with respect to the corresponding action variable, 
hence they can be recognized as the radial {ujr) and angular {uu,^) frequencies of the classical 
orbits respectively. Thus 

e(n^, I') - e{nr, I) ~ (n^ - nr)oUr + {I' - O'^v' > (2-2) 

(we use units such that h = c = 1). This combination of frequencies is exactly that appear- 
ing in the denominator of the zero-order response function obtained in [2] from the solution 
of the linearized Vlasov equation with fixed-surface boundary conditions. Moreover, from a 
comparison of the zero-order Vlasov response function with the analogous quantum propa- 
gator [8], it can be seen that in the Vlasov propagator the quantum matrix elements of the 
excitation operator are replaced by appropriate Fourier coefficients that can be evaluated 
as integrals along the classical trajectories. Since here we are interested in the octupole 
response, we report the propagator describing the response of a nucleus to an external field 
of the type 

g(r) = r-V3M(f) (2.3) 
in the single-particle approximation. This propagator is given by [5] 



^Ls « = ^ E E C^N / dxx^ SnN{x)^—- f-. 2.4 
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Instead of the frequency a;, as independent variable we have used the dimensionless 
quantity s = u / {vp/ R), as a, consequence the A-dependence of this propagator is factorized 
as AR^ oc A^. The eigenfrequencies (2.2) are accordingly replaced by the functions (for a 
square-well potential) 

SnN{x) = ^. 2.5 

X 

The variable x is related to the classical angular momentum A of a nucleon. The relation 
is x = sin a, where a is the angle spanned by the radial vector when the particle moves from 
the inner to the outer turning point. Clearly, for a square- well potential one has cos a — X/X, 
where A is the maximum particle angular momentum A = ppR, is the Fermi momentum, 
while Vp and ep are the corresponding velocity and energy. The sum over discrete angular 
momentum values of the quantum propagator is replaced here by the integration over x, 
thus the particle angular momentum is treated as a continuous variable. 

The quantities Csn in Eq. (2.4) are classical limits of the Clebsh-Gordan coefficients 
coming from the angular integration of the quantum matrix elements [8]. Their explicit 
value is Cs±i = ^ and Cs±s = In principle the integer N takes values between —L and 
L, however only the coefficients CpN where has the same parity as L are nonvanishing. 
The coefficients Qnl/i^) appearing in the numerator of Eq. (2.4) have been defined in Ref. 
[2], they are essentially the classical limit of the radial matrix elements of the octupole 
operator and can be evaluated explicitly as 



For = ±1 this expression coincides with that appearing in the (compression) dipole 
response (cf. Eq.(A.5) of Ref. [5]), however for the octupole response we also need terms with 
N — ±3. These new modes have an interesting property since the associated eigenfrequencies 
SnN{x) can vanish in the interval < a{x) < | [the equation {nn + Na) = has a solution 
for a = |, corresponding to closed triangular orbits]. In Ref. [1] it has been pointed out 
that the vanishing of this eigenfrequency might give rise to a possible instability against 
octupole-type deformations in nuclei heavier than ^°^Pb. Although at first sight it might 
seem that the coefficients (2.6) would diverge when SnN{x) — > 0, it is actually possible to 
check that 

liIn_gg±3(a;) = -li?^ (2.7) 

The very fact that this limit is finite is important for our discussion about the role of 
triangular nucleon orbits. The linear response function in the single-particle approximation 
is well-behaved because 

lim(Im7^°(s)) =0, (2.8) 

and 

limfRe7^°(s)) = ^A^ , (2.9) 
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a result that we shall need later. Since the zero-order octupole strength function is propor- 
tional to the imaginary part of the response function (2.4), we can see that the contribution 
of triangular trajectories to the octupole strength function does not give any pathology, at 
least at the single-particle level. Of course instabilities could still arise in the collective re- 
sponse, so we study also this response by using the same model, however, before doing this 
we look more in detail to the single-particle octupole strength given by the Vlasov theory. 
The zero-order octupole strength function Sl=3{E), given by {E — huj) 

5i=3(^) = --Im7e°(E) , (2.10) 

TT 

is shown in Fig.l. It can be seen that the single-particle octupole strength is concen- 
trated in two regions around 8 and 24 MeV. As pointed out already in [2], in this re- 
spect our semiclassical response is strikingly similar to the quantum response, which is 
concentrated in the ITiuj and S^tu regions. The modes that contribute most are those 
with (n,iV) = (0,1), (1,-1), (-1,3), (2,-3), for the lliuj region and with (n, iV) = 
(0,3), (1,1), (2,-1), (3,-3), for the "ihu region. This concentration of strength is quite 
remarkable because our static distribution, which is taken to be of the Thomas-Fermi type 
/o(r, p) oc 9{ep — /io(r,p)), does not include any shell effect, however, because of the close 
connection between shell structure and classical trajectories expressed by Eqs. (2.2) and 
(2.6), we still obtain a stregth distribution that is very similar to the one usually attributed 
to shell effects. Clearly the integers n and N correspond to the difference of radial and 
angular momentum quantum numbers in Eq. (2.2). 



III. COLLECTIVE RESPONSE (FIXED SURFACE) 

The zero-order response function (2.4) gives only a first approximation to the nuclear 
response. When the residual interaction between nucleons is taken into account, a collective 
response function can be obtained by solving the Vlasov equation with appropriate boundary 
conditions [2,3]. If the interaction is assumed to be of the octupole-octupole type, 

y(ri,r2) = K.rlrl Y.Y;^{v,)Y^M{i2) (3.1) 

M 

the collective fixed-surface octupole response function is given by [9] 

The parameter specifies the strength of the residual interaction, its value can be estimated 
in a self-consistent way, giving ( [1], p. 557), 

KBM = -f ^ ^ -1- lO-^MeV/fm^ (3.3) 

for the octupole case. The parameter ujq is given by ujq ^ AlA'^MeV . Since this es- 
timate is based on a harmonic oscillator mean field and we are assuming a square-well 
potential instead, we expect some differences. Hence we shall determine the parameter 
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phenomenologically, by requiring that the peak of the high-energy octupole resonance agrees 
with experiment. This requirement impUes K3 ^ 2kbm, which is in agreement with the pre- 
scription obtained in the quadrupole case [6]. In Fig.2 we report the coUective octupole 
strength function given by Eq. (3.2), with this value of K3. We can clearly see the effects 
of collectivity that result in a shift and concentration of the strength into two sharp peaks 
around 20 Mcv and 6-7 Mev. The experimentally observed [10] concentration of isoscalar 
octupole strength in the two regions usually denoted by HEOR (high energy octupole reso- 
nance) and LEOR (low energy octupole resonance) is qualitatively reproduced, however the 
considerable strength experimentally observed at lower energy (low-lying collective states) 
is absent from our fixed-surface response function. Like for the quadrupole response [6], we 
need to consider a different solution of the linearized Vlasov equation (in which the nuclear 
surface is allowed to vibrate [3]) in order to account for this feature of the response. 

The poles of the collective response function are determined by the vanishing of the de- 
nominator in Eq. (3.2). A solution of this equation at zero or purely imaginary frequency 
could be interpreted as an instability of the spherical shape against an octupole-type defor- 
mation of the ground state, however we can see that this could happen only for interactions 
stronger than 

i^scr = ~ -3.14 lO-^MeV/fm^ . (3.4) 

Our value of K3 is smaller (in absolute value) than this, so in the present model the spherical 
shape is stable. This is vahd for any value of A, as long as we take = 2k,bm: since the 
product Ksm'^3(0) is A-independent. 

The triangular orbits do not enter directly into the discussion of this collective response, 
but the fact that they do not generate pathologies at the level of the single-particle response 
7^.3(5) simphfies the discussion of the stabihty in the collective response also. 



IV. COLLECTIVE RESPONSE (MOVING SURFACE) 



The collective fixed-surface response function (3.2) docs reproduce two important fea- 
tures of the octupole response: the HEOR and the LEOR, however it misses the experimen- 
tally observed low-lying collective 3^ states. In order to account for these low-lying states, 
a different solution of the linearized Vlasov equation has been proposed in Ref. [3]. In that 
approach the nuclear surface is allowed to vibrate according to the usual liquid-drop model 
expression 

R{e, ip,t) = R + Y, SRLM{t)YLM{e, if) . (4.1) 

LM 

A self-consistency condition involving the surface tension is then used to determine the time- 
dependence of the additional collective variables 5Rlm{^)- Then, always for a separable 
residual interaction of the kind (3.1), the collective fixed-surface response function (3.2) is 
replaced by the following moving-surface response function [6] 

^3(s)=7^3(s)+53(s), (4.2) 
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with 7^3(5) still given by Eq. (3.2), while 5,3(5) gives the moving-surface contribution. For 
the separable interaction (3.1) the function S^ls) can be evaluated explicitly as [6] 

P., [xlis) + K,goR'nl{s)r .... 

with C3 = lOaR^ + {C3)coui (c" ~ IMeVfm"^ is the surface tension parameter obtained 
from the mass formula, (6*3)20^; gives the Coulomb contribution to the restoring force) and 
Qo — A/^R^ the equilibrium density. 

The functions Xsi^) ^^'^ Xai^) defined as 

XIW = M 2 r (4.4, 



and 

their structure is similar to that of the zero-order propagator (2.4). 

In Fig. 3 we report the octupole strength function given by the moving-surface response 
function (4.2) and compare it to the collective fixed-surface response given by Eq. (3.2). The 
most relevant change induced by the moving surface is the large double hump appearing at 
low energy. This feature is in qualitative agreement both with experiment [10] and with the 
result of RPA-type calculations [11,12], moreover is rather similar to that found in Ref. [13] 
within the same model, but with different excitation operator and residual interaction. We 
interpret this low-energy double hump as a superposition of surface vibrations and LEOR. 
In order to explain why we do not obtain one or more sharp 3' states at low energy, we 
have to analyze our moving-surface response function in some detail. Since the explicit 
expression (4.3) looks rather involved, we consider the limit of non-interacting nucleons. If 
we let K,3 — > 0, the function (4.3) becomes 

Slis) = -i?^ , (4.6) 

and the full moving-surface response function (4.2) 

nlis) = n',{s) + slis) . (4.7) 

In Fig. 4 we show the octupole response function evaluated both for ^3 = 2kbm (solid 
curve) and for K3 = (dashed curve). We can see from this figure that, while the residual 
interaction changes drastically the response in the giant resonance region, at low energy the 
octupole response is affected only very slightly by the interaction (3.1). Thus the low-energy 
octupole response can be analyzed by using the simpler formula (4.6), rather than Eq. (4.3). 

The eigenfrequencies of low-energy collective modes are approximately determined by the 
vanishing of the denominator in Eq. (4.6), moreover at low frequency the function Xai^) 
can be expanded as [4] 
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(4.8) 



implying that the parameters 6R3M{t), that describe octupole surface vibrations in Eq. 
(4.1), approximately satisfy an equation of motion of the damped oscillator kind: 

DsSR^uit) + 73SR3M{t) + CsSRsMit) = . (4.9) 

It has already been pointed out in [4] that for A = 208 the numerical values of parameters 
in this equation are such that this oscillator is actually overdamped, moreover, always in [4] , 
the coefficients have been evaluated analyticallty in the low-frequency limit, giving (for 
a generic L) 

= ^wf '^ 2L+1 ^ N ' ^^^^2' 2^ I]cosanivsin2aniv0(- - ^niv) , (4.10) 

with 7^/ = ^qqPpR'^ and anN = Ji^^- The angles anN are related to the nucleon trajectories 
as discussed in Sect. II. In our case the coefficient ^l=z gets a contribution only from the 
term with n—1 and = 3. thus we see that only nucleons moving along closed triangular 
trajectories can contribute to the damping of octupole surface vibrations. In order to check 
that indeed the closed triangular orbits are the main source of Landau damping in the low- 
energy octupole response, we study also the response that is obtained when the contribution 
of these orbits is excluded. This exclusion can be made simply by avoiding a small interval 
of length 2 A about the value a; = ^ in the x- integration in Eqs. (2.4), (4.4) and(4.5). In 
Fig. 5 we show the response obtained with a value of A = 0.02. The shape of the low-energy 
hunp changes dramatically because of the lack of damping due to the missing triangular 
trajectories. A very sharp peak is now developed at low energy. The position of this peak 
is influenced by the value of the mass parameter D^, omitting the triangular orbits changes 
also the numerical value of this parameter. Once again we see that the closed triangular 
orbits, rather than generating a shape instability in the octupole channel, are the main 
source of Landau damping in this response function. Of course quantum effects can play an 
essential role at this level since in quantum mechanics the angular momentum is quantized 
and there could be no value of angular momentum corresponding to triangular trajectories. 

We would also like to point out that our moving-surface response function can display 
a shape instability in the octupole channel, which can arise if the restoring force parameter 
C3 vanishes. This could happen if the repulsive Coulomb term {C3)coui exactly balances 
the attractive surface-tension part of the restoring force. In this case our moving-surface 
response function developes a pole at the origin, as shown in Fig. 6. For odd multipolarities 
the instability condition in our model is the same as for the liquid-drop model [4] . 

Before concluding, in Fig. 7 we compare our moving surface response function with the 
response function of the overdamped oscillator (4.9). This figure supports our interpretation 
of the lower-energy hump as due to overdamped surface vibrations. 

V. CONCLUSIONS 

We have analyzed the octupole response function of a hypothetical heavy nucleus con- 
taining A — 208 nucleons in a square-well potential by employing a semiclassical model 
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that relates the response to features of the classical trajectories of nucleons. Triangular tra- 
jectories are particularly relevant for the octupole response because, at the single-particle 
level, a simple combination of their characteristic frequencies, appearing in this response 
function, can vanish. This fact had already been noticed long ago and some speculations 
had been made on the possible connection between these classical orbits and the onset of 
a shape instability in heavy spherical nuclei against octupole-type deformations [1]. Our 
detailed calculations performed in a semiclassical approach show that the vanishing of the 
eigenfrequencies associated with triangular trajectories has no consequences on the stability 
of the spherical shape in heavy nuclei, rather, triangular orbits are essential in providing a 
damping of low-energy octupole excitations. 
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FIGURES 



Fig.l Semiclassical octupole strength function analogous to quantum single-particle 
strength function. Note the strength concentration in the Ihcu and 3hu! regions. Calcu- 
lations are for A — 208 nucleons in a square well potential of radius R — 1.2A3 fm. 

Fig. 2 Collective strength function (solid curve) evaluated for fixed-surface boundary con- 
ditions. The residual interaction parameter is % = —2. 10~^MeV/fm^. The dashed curve is 
the strength function of Fig.l. 

Fig. 3 The solid curve shows the octupole strength function evaluated by using moving- 
surface boundary conditions in the Vlasov equation, the dashed curve instead, corresponds 
to fixed-surface boundary conditions. 

Fig. 4 Moving-surface strength function for = —2. 10~^MeV/fm^ (solid) and for = 
(dashed). 

Fig. 5 Moving-surface strength function with (solid) and without (dashed) contribution 
of closed triangular orbits. 

Fig. 6 Collective moving-surface strength function (solid) dispaying a divergence at van- 
ishing excitation energy when C3 ^ (dashed). This divergence is interpreted as a shape 
instability. 

Fig. 7 Collective moving-surface strength function (solid) compared to overdamped oscil- 
lator strength function (dashed) with appropriate parameters. The two functions practically 
coincide for < 2 MeV. 
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